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Abstract. We investigate random walks on a lattice with imperfect traps. In one dimension, we pertur- 
batively compute the survival probability by reducing the problem to a particle diffusing on a closed ring 
containing just one single trap. Numerical simulations reveal this solution, which is exact in the limit of 
perfect traps, to be remarkably robust with respect to a significant lowering of the trapping probability. 
We demonstrate that for randomly distributed traps, the long-time asymptotics of our result recovers the 
known stretched exponential decay. We also study an anisotropic three-dimensional version of our model, 
where for sufficiently large transverse diffusion the system is described by the mean-field kinetics. We 
discuss possible applications of some of our findings to the decay of excitons in semiconducting organic 
polymer materials, and emphasize the crucial influence of the spatial trap distribution on the kinetics. 

PACS. 02.50.Ey Stochastic processes. - 05.40.-a Fluctuation phenomena, random processes, noise, and 
Brownian motion. 05.60.-k Transport processes. 



1 Introduction 

Random walks in disordered media H,0| form a class of 
statistical processes that have been widely employed for 
the description of an impressive number of physical, chem- 
ical and biological phenomena. A particularly interesting 
subclass is constituted by Brownian motion in the pres- 
ence of quenched disorder in the medium. Depending on 
the nature of this disorder, significant deviations from the 
results known for pure random walks are observed. This 
is true for microscopic quantities (e.g. first passage times, 
average distance from origin at time i, etc.) as well as 
for macroscopic properties, such as the asymptotic time 
evolution of the overall population of random walkers un- 
dergoing site-correlated annihilation [Q,^. 

In the present paper, we address the intermediate and 
long-time kinetics of random walkers decaying through 
two simultaneous channels: (i) spontaneous ('radioactive') 
decay, with no spatial or temporal correlations involved, 
and (ii) a decay process induced by the presence of non- 
moving (quenched) imperfect traps in the medium. 

In the existing literature, the term 'trap' has been used 
to describe different types of disorder: In certain instances, 
the concept of a trap has embodied potential wells in the 
medium, acting on a walker through delaying its motion 
(the 'valley model') 0]. In the present context, we re- 
fer to traps simply as sites where particles may undergo 



spontaneous annihilation with a certain (fixed) probabil- 
ity < g < 1. 

Random walks with annihilation by trapping have al- 
ready been the object of considerable attention during the 
past two decades |,|,|,0,|,|,|l^,0|l|, to name but a 
few which are of relevance for this work. These works rep- 
resent a fundamental and necessary step towards an un- 
derstanding of transport processes in partially absorbing 
media. In biology, this applies, e.g., to the scattering of 
laser light in tissues. In chemistry, microscopic theories 
of chemical reaction kinetics necessitate to access quan- 
tities such as the probability distribution of the nearest- 
neighbor distance of a diffusing molecule to a (static or 
moving) trap, representing a reaction center. 

In solid-state physics, the study of a whole set of dif- 
ferent phenomena in disordered systems relies heavily on 
the understanding of random walks in the presence of 
partially absorbing traps: Electron-hole recombination in 
amorphous solids, and chemical binding by impurities of 
interstitial hydrogen atoms in metals represent two com- 
mon examples. More indirectly, this problem is also re- 
lated to self-attracting polymers, and to the density of 
states of binary disordered systems, see Ref. [Q. 

We draw part of our motivation for this present study 
from a class of organic semiconducting materials which 
has been the object of considerable attention recently [p^ . 
Due to the nature of their opto-electronic properties, thin 
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films of organic molecules working as semiconducting dio- 
des ||l4| are already employed in high-performance elec- 
troluminescent displays |15| . For the initial experimental 
studies of excited states in these materials, aiming at un- 
derstanding their fundamental properties, light is used for 
their generation — in most cases the elementary excita- 
tions take the form of a bound electron-hole pair |16t . 
These excitations behave as pseudo-particles (excitons) 
which diffuse and interact in varying and complex ways 
with the medium p7| , and eventually annihilate by emit- 
ting photons or phonons. The latter, non-radiative anni- 
hilation process is known to be mediated by chemical de- 
fects |l|. 

Experimental studies suggest that the physics of re- 
laxation at work in such excitable media depends in a 
rather tight and non-trivial way on a variety of partially 
controllable parameters. Among them, we quote the chem- 
ical nature of the material used, the way it has been pre- 
pared and/or altered by techniques like photo-oxidation, 
and the configurational structure at the microscopic level 
as well as interactions between individual polymer chains. 
These mutual interactions arc markedly different in solid 
or liquid solutions [[l9| , and amorphous or polycrystalline 
films, respectively. The detailed relaxation mechanisms, 
and their dependence on sample properties are as yet 
incompletely understood, as far as their relation to the 
elementary and fundamental processes involved are con- 
cerned. We believe that the present, somewhat unsatis- 
factory situation calls for an effort in investigating, e.g., 
through simple kinetic statistical models, the role played 
by these different parameters. 

In addition, the effectively one-dimensional (or at least 
strongly directed, anisotropic) nature of cxciton transport 
in organic semiconductors renders simple mean-field the- 
ory approaches inadequate. Rather, one expects marked 
fluctuation and correlation effects. For the model studied 
here, these in turn should strongly depend on the spatial 
defect distribution. Thus, experiments probing the exci- 
ton kinetics in semiconducting polymers promise to pro- 
vide an excellent testing ground for a wide variety of non- 
equilibrium statistical models. Indeed, the dynamics of 
laser-induced excitons in N(CIl3)4MnCl3 (TMMC) poly- 
mer chains appears to be the best experimental realization 
to date for diffusion-hmited fusion processes A + A ^ A, 
where (in an intermediate time window) the predicted 
power law for the particle density n ~ t~^/^ has actually 
been observed unambiguously [Q. 

The present paper is organized as follows: In Sec. |^ we 
provide a precise definition of our model of a random walk 
with imperfect traps (inducing non-radiative exciton re- 
combination) and spontaneous annihilation (or radiative 
recombination). Our main interest being the intermediate 
and long-time decay of the overall population of particles, 
the following section presents a calculation of the tempo- 
ral behavior of this quantity in a one-dimensional setup. It 
is known that the low dimensionality has a strong impact 
on the late-time asymptotics of the decay rate [|[|3[^ . Our 
analytical method is based on a scheme that we name the 
'decoupled-ring' approximation. Similar approaches have 



been used in the literature ilO,p[ . We demonstrate that the 
late-time resummation of our analytical result reproduces 
the known asymptot ic, very slow stretched-exponential 
decay |@,||,||]. Section 4J investigates, by means of Monte 
Carlo simulations, the range of validity for this analytic 
result when the trapping probability departs from one, 
the value for which th e de coupled-ring approximation be- 
comes exact. In Sec. 4.2, we examine, by way of addi- 



tional simulations, how the results presented resist to a 
relaxation of the one-dimensional constraint, by allowing 
particles to slightly diffuse also in the two transverse di- 
rections. Finally, in Sec. |^ we discuss the possibility of 
applying such a model to account for experimental mea- 
surements of exciton decay in organic semiconductors, and 
comment on the crucial role of the spatial trap distribu- 
tion for the long-time kinetics. 



2 The model 

We consider a one-dimensional, regular, infinite lattice C 
with linear lattice spacing a. A fraction c of the lattice sites 
arc tagged, thereby denoting that they possess a special 
property. These particular sites will hereafter be referred 
to as 'traps', indicating their specific role in the dynamics 
we are about to introduce. 

In general, the spatial trap distribution can be arbi- 
trary. We shall single out two generic, contrasting situa- 
tions: 

1. The trap distribution is random, i.e., the distance be- 
tween two consecutive traps along the chain is repre- 
sented by a stochastic variable with poissonian distri- 
bution. 

2. The location of the traps is spatially correlated. This 
case typically covers the situation of n-modal intertrap 
spacings, or, more generically, situations where some 
kind of regularity (order) can be identified in the way 
traps are distributed over the system. 

Both situations are encountered in experimental samples: 
If properly synthesized, semiconducting polymer chains in 
solution may be considered to be essentially trap-free, ex- 
cept for their extremities perhaps, thus falling into the sec- 
ond category. Photo-oxidation of the same material leads 
to the creation of defects at random locations along the 
chain, and the first of the two above cases then applies. 

We now introduce entities, representing (pseudo- or 
quasi-)particles, initially randomly distributed over the 
lattice sites, with density < po < 1- A discrete dynamics 
is implemented, by allowing particles to perform symmet- 
ric random walks: A particle located at site k at time t can 
subsequently be found either at (fc — 1, i -I- r) with proba- 
bility 1/2, or at (k + l,t + r) with the same probability, 
where r denotes the physical timestep associated to the 
iteration of the dynamics. 

At each timestep, every diffusing particle may annihi- 
late through two different decay channels: 

— Spontaneous decay, occuring with fixed probability < 
p < 1. For excitons in organic semiconductors, this 
channel represents radiative recombination. 
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— Stochastic capture by a trap: If the diffusing particle 
enters a tagged site, it annihilates with a given prob- 
ability < g < 1. For our later application to exciton 
kinetics, we specify this trap-mediated decay process 
to be non-radiative, i.e., not emitting light at visible 
wavelengths. 

In both cases, the particle is simply removed from the 
system. The situation with q = I will be referred to as 
perfect traps. 

As both decay process act concurrently without di- 
rectly influencing each other, their combined effect ap- 
pears as the product of two independent factors. The spon- 
taneous decay being trivially described by a simple expo- 
nential relaxation, the real difficulty here rests in the con- 
tribution to the kinetics related to the traps, where the dy- 
namical build-up of spatio-temporal correlations renders 
any mean-field type of approach invalid in (sufficiently) 
low dimensions. As we shall see, the spatial arrangement 
of the traps may crucially determine the long-time asymp- 
totics of the particle decay in this situation. 



3 Asymptotic kinetics in the decoupled-ring 
approximation 




Fig. 1. Schematic iUustration of the ring -decoupling procedure. 
Black dots represent regular lattice sites, white ones those with 
a trap. The infinite, one-dimensional chain is split into seg- 
ments (delimited by the location of the traps), each of which 
is then closed onto itself to form a ring. The entire chain is 
thereby mapped onto an infinite set of rings of varying lengths. 
This approximation becomes exact in the limit of perfect traps. 



3.1 Motivations 



3.2 Single-ring solution 



The model of perfect, uncorrelated traps in one dimension 
has been solved by Anlauf ||] by means of the span dis- 
tribution function for one-dimensional random walks. The 
model of imperfect trapping has been treated in the liter- 
ature by various means. Weiss and Havlin Q studied the 
system in one dimension by introducing a modified model 
where particles are either destroyed or reflected at a trap, 
thus decoupling the line segments separated by a trap. 
This calculation was redone and some errors were cor- 
rected in g . Additionally, the coupled system was mapped 
to a harmonic chain with random masses and the authors 
were able to obtain the long-time behaviour of the parti- 
cle decay in one dimension. We chose a similar approach 
to ||10| by decoupling line segments by cutting the infinite 
line at each trap and bending each segment onto itself to 
form a ring (see Fig. m). The advantage of this method 
over the approach of mapping to a harmonic system with 
random masses is that it is easy to obtain results for cor- 
related traps: since each decoupled ring can be treated 
seperately, the final result is simply obtained by averag- 
ing over the ring lengths with the appropriate probability 
distribution. The disadvantage is that it contains uncon- 
trolled approximations since particles can never leave their 
ring and wander off into other regions as they could in 
the original model. However, computer simulations show 
that these approximations are surprisingly good, in par- 
ticular for not too low trapping probabilities. For perfect 
traps, the decoupled model even becomes exact since the 
traps prevent passage between segments even in the orig- 
inal model. 



Here, we briefly sketch the solution of the model on a sin- 
gle ring of length n. The discrete time dynamics takes the 
form of a master equation for the particle density distribu- 
tion vector p{t) = {pi{t), . . . , Pn{t)), with the trap located 
at site 1 without loss of generality. 
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For late times, the overall decay rate of the concentra- 
tion is 'slaved' to the largest eigenvalue a„ of the above 
matrix, and its corresponding eigenvector a. T he ca lcu- 
lation is straightforward and is shown in App. |A.l| ; the 
result is 



a„ = cos (pn , where 
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for the eigenvector. 

Starting from an initial distribution pk{t = 0), the 
long-time solution is 



reads: 



, . t,n^oo {a,p{Q)) ^t\nc 

[a, a) 



(6) 



where (•, •) denotes the scalar product. The scalar prod- 
ucts can be worked out for a homogeneous initial distri- 
bution, Pfc(O) = poi to leading order in 1/n, and the result 
for the mean density on the ring p is 



1 " 

Tt ^ ^ 



k=l 






(7) 



Note that this result is better than might be expected at 
first sight because the overlap of the eigenvector a with 
the i nitia l homogeneous configuration is very large (see 
App. Q, 



(a,p(0)) 2^/2 



I«IIP(0)| 



+ O{l/n)K0.9. 



(8) 



which guarantees that this result is good even at interme- 
diate times. 



3.3 Full-chain solution 

Given an arbitrary distribution of traps, one can now write 
down the solution for the full chain (in the decoupled-ring 
approximation) by simply superposing the single-ring re- 
sults, properly weighted with the relative abundance for 
the occurrence of rings of length n in the chain. For exam- 
ple, if there are well-defined trap spacings such that only 
a finite number of certain values of n^, i = 1, . . . , /c is pos- 
sible, with relative abundance P{ni), then the final result 
is simply a sum of k terms of the form IMl each multi- 
plied with P{ni). As is evident from Eq. (R), the terms 
for large n will quickly dominate, and the maximum value 
'T-max will govern the long-time limit. 

For a random (poissonian) initial distribution of the 
traps with concentration < c < 1, all values of n are 
possible, the probability of finding a ring of size n being 
given by 



P(n)=nc2(l-c)"-i . 



This allows us to finally write 
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For i — > oo, the late-time asymptotics of the infinite 
sum ( llOl ) can be extracted (sec the appendix). The result 
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We recover the known stretched exponential long-time 
decay [§|,|[|^ p(i) ~ exp(— const, i^/'^) with the correct 
values for both the exponent and its prefactor, and in- 
cluding a g-dependent enhancement factor in agreement 
with B (higher order correction terms have been omitted 
here for brevity). With a random (Poissonian) distribution 
of traps, one therefore expects a much slower asymptotic 
time decay than for a typical correlated distribution. Cer- 
tainly, for any finite number of allowed values of n, the 
asymptotic temporal decay will be a simple exponential. 
As another example, a Gaussian distribution of the lengths 
of trap-free regions (around some mean value) leads to a 
stretched exponential with exponent 1/2 (by a calcula- 
tion analogous to the one above). Thus, in addition to the 
trap concentration and the induced decay probability g, 
the spatial arrangement of the traps is of crucial impor- 
tance. 

Instead of using the asymptotic result Eq. (|ll]) we 
prefer to use Eq. ( |lO| ) (suitably truncating the sum) for 
our comparisons with simulations and experiments below 
since it applies not only to the long times but also to in- 
termediate times. This is essential as it is known that the 
crossover to the long-time behaviour Eq. (O) may set in 
so late to be unreachable in practice (see e.g. [Q and 
references therein, |7|,[12[). 



4 Monte Carlo simulations 

4.1 Strictly one-dimensional system 

We have investigated the range of validity of the result 
(nOf) emerging from our decoupled-ring approximation by 
means of Monte Carlo (MC) simulations. The parameter 
of interest is of course g, the trapping probability. For q 
equal (or very close) to 1, Eq. ( |lO|) should perform very 
well. This agreement may be expected to degrade when q 
is lowered significantly below unity, since this implies an 
increasing coupling between adjacent rings that has been 
completely neglected in our treatment of the problem. 

Figure shows a set of plots comparing results of MC 
simulations to the approximate analytic formula dlQ), for 
different values of the trapping probability q. The spa- 
tial trap distribution is random (poissonian). For clar- 
ity, we have set the spontaneous decay probability p to 
zero here. No crossover may anyway be expected between 
the exponential and stretched exponential components of 
the dynamics: Due to the structure of Eq. ( pi] ) , the sub- 
exponential behavior is screened out by the spontaneous 
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Fig. 2. Monte Carlo simulations for the one-dimensional im- 
perfect trapping process. The dimensionless quantity on the 
vertical axis is the amount of random walkers present in the 
system, normalized by the initial total population. It can also 
be interpreted in a statistical sense, as the survival probability 
of a single walker. Crosses represent the numerical result, aver- 
aged over ~ 5000 independent runs. The solid lines represent 
the analytical result of the ring-decoupling scheme, as given by 
Eq. (|lO|). The following choice of parameters was made: system 
size=2 '^ — 16384 sites; p (spontaneous decay probability) — 0; 
c (trap concentration) — 0.2; diffusion coefficient — 0.5 [a^ /t]; 
po (initial particle density) = 0.6. 



exponential decay for all real and positive values of the 
time t. 

The first (upper left) graph for q = 1 serves as a ref- 
erence for a visual evaluation of the uncertainty in our 
numerical data. Looking at the other graphs in the fig- 
ure, one observes that for the chosen trap concentration 
(c = 0.2), the analytical result performs remarkably well 
down to q = 0.4 (!). Note that the finiteness of the systems 
simulated implies a temporal horizon: Beyond this limit, 
both strong depletion in particles and inaccurate statistics 
for large, trap-free regions (which play an increasingly im- 
portant role for late times) conspire to render the simula- 
tion meaningless. Our approximation being an asymptotic 
theory, there is however no doubt that the matching be- 
tween Eq. dlOl ) and the numerical data is at least as good 
for the unexplored, late time domain as it appears in the 
'early-time' regime depicted in Fig. y. 

Figure S illustrates how the decoupled-ring approxi- 
mation performs upon varying the trap concentration at 
fixed trapping probability. 



4.2 Anisotropic three-dimensional system 

Semiconducting polymer chains can be found in various 
configurations: When in solution, they can be considered 
as totally isolated from each other, hence representing a 
set of one-dimensional segments on which excitons propa- 
gate and interact. They can, however, also be organized in 
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Fig. 3. Monte Carlo simulations for the one-dimensional im- 
perfect trapping process with fixed trapping probability and 
varying traps concentration. The quantities on the axis and 
the meaning of the symbols chosen are the same as for fig- 
ure 0. Results have been averaged over ~ 5000 independent 
runs. The following choice of parameters was made: system size 
= 2^* — 16384 sites; p (spontaneous decay probability) — 0; q 
(trapping probability) = 0.4; diffusion coefflcient = 0.5 [a /r]; 
Po (initial particle density) = 0.6. 



a more cohesive manner, namely in the form of a network 
characterized by a certain degree of coupling between the 
polymer chains, or even in highly ordered structures such 
as polycrystalline films. This interaction between chains 
is evident both in the chain-to-chain transfer of charges 
when these networks are used in light emitting diodes, as 
well as in the chain-to-chain transfer of excitons [g2|. In 
these cases, the one-dimensional character of the exciton 
dynamics is likely to be broken, the 'particles' being able 
to 'cross-diffuse' from one chain to a neighboring one. 

In this context, an interesting question that naturally 
arises is the following: To which extent does the picture 
related to the ideally one-dimensional case break down, 
once a small degree of cross-diffusion is allowed in a set 
of one-dimensional chains forming a crystal- like structure? 
For we know that in an isotropic three-dimensional system 
with uncorrelated traps, the long-time behaviour of the 
system is given by a stretched exponential of the form 
p{t) - exp(-const.t3/5) [||]. 

Our one-dimensional model is extended to three di- 
mensions by allowing particles to hop in all three spa- 
tial directions, but with different rates. In particular, we 
single out one direction as the direction of the polymer 
chains and allow hopping along this direction as in the one- 
dimensional model with a diffusion constant Di. Trans- 
verse hopping is suppressed by choosing a transverse dif- 
fusion constant Dt < Di. Thus, a particle hops along the 
polymer direction with probability Pi = {2Dt/Di + 1)^^ 
and perpendicular to it with probability 1 — Pi at each 
time step. 

Figure H displays the result of MC simulations of our 
model in three dimensions, where the ratio of the trans- 
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Fig. 4. Monte Carlo simulations of a three-dimensional imper- 
fect trapping process competing with spontaneous decay (see 
the caption of Fig. H for comments on the axis and units). 
Symbols represent the numerical results for different values of 
the transverse diffusion coefflcient Dt, as indicated in the leg- 
end. They have been averaged over ca. 100 runs. The solid line 
shows a fit with a stretched exponential exp(— ai^' ^) as appro- 
priate for a truly three-dimensional system, showing that at 
Dt/ Di — 0.002 and c = 0.2 the system already behaves like a 
three-dimensional one. The following choice of parameters was 
made: system size = 2^^ x 2^ x 2^^ sites; p = 0.001; c = 0.2; 
q — 0.5; longitudinal diffusion coefficient Di = 0.5 [a^ /t]; 
po = 0.6. 



verse and longitudinal diffusion rates Dt/ Di has been cho- 
sen non-zero, but kept small. The virtue of this figure is 
to undoubtedly show the extent to which the topology 
affects the dynamics in this process: Even a very small 
relaxation of the onc-dimcnsional constraint results in a 
drastic effect. This is easily understood if one takes into 
consideration the two following important aspects of the 
dynamics: 

— The increasing role played by trap- free regions of larger 
and larger size as time increases. 

— The fact that a particle jumping transversely onto a 
neighboring chain has, statistically, large chances to 
penetrate a small inter-trap segment, and therefore de- 
cay rapidly. This is of course true only under the hy- 
pothesis that the trap distributions for all chains are 
totally uncorrelated with respect to each other. 

On the other hand, we see that the decoupled-ring ap- 
proximation, which correctly captures only the one-dimen- 
sional topology, numerically remains a fairly adequate de- 
scription for small anisotropy ratios. 

In Fig. H, we examine how the trap concentration af- 
fects the amplitude of the effect of the relaxation of the 
one-dimensional constraint on the population decay. As 
shown, a higher trap concentration helps to screen out the 
emergence of the transverse channel for diffusion. This is 
an important observation when one tries to use our model 
to interpret experimental data, as we will see in the forth- 
coming section. 
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Fig. 5. Effect of the trap concentration on the relaxation of 
the one-dimensional constraint. The same simulations as in 
Fig. M were performed for a higher trap concentration c = 0.7 
(averaged over ~ 100 runs here, too), and both sets are now 
displayed together. Solid lines correspond to c = 0.2, and dot- 
ted ones to c = 0.7. For each of the three cases with Dt 7^ 0, 
the plot shows the relative deviation from the Dt = case. 
Circles, squares, and triangles indicate respectively the simu- 
lations with values Dt/Di = 0.002,0.02, and 0.2. 

5 Applicability to the decay of excitons in 
quasi-one-dimensional organic 
semiconductors 

As previously mentioned, our model may be considered a 
potential candidate for a 'minimal' description of the de- 
cay dynamics of excitons in organic semiconductors. Con- 
frontation with real data coming from experimental mea- 
surements of such a decay is therefore instructive. The fig- 
ures below display our attempts to fit such experimental 
data with our analytical result (|lO|), supplemented with a 
spontaneous decay factor e^^*. 

The data presented here, orginally reported in Rcf. [£3 , 
measure (in frequency space) the response of ladder-type 
Poly(Para-Phenylene) (LPPP) to excitation by a modu- 
lated light beam. The signal represents the Fourier trans- 
form of the spectrally integrated photoluminescent emis- 
sion of the material. This is given by the activity, i.e., the 
number of radiative recombinations of excited electronic 
states (excitons) per time interval. Thus, while both spon- 
taneous and trap-induced decay channels reduce the exci- 
ton population, only the activity due to the spontaneous 
radiative recombination is monitored by the data. Three 
different types of samples, all made from the same ma- 
terial, were analyzed: The polymer in solution, a pristine 
polymer film, and the same film after it had been sub- 
jected to photooxidization. In the following, we report our 
comparison of theoretical and experimental data for each 



Figures g-jg show the experimental data as real and 
imaginary parts. The solid lines are simultaneous fits (per- 
formed with a Levenberg-Marquardt algorithm) to both 
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Fig. 6. Fourier transform of the response function as measured 
(dots) and a fit (solid lines) with the Fourier transform of the 
activity for an exponential time decay for the polymer in solu- 
tion. Fit parameters are the spontaneous decay rate p and an 
(irrelevant) overall amplitude A, given in arbitrary units. 
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Fig. 8. Fourier transform of the response function as measured 
(dots) and a fit (solid lines) with the Fourier transform of the 
activity corresponding to the time decay according to En. ( |l0[ ) 
for the non-oxidized film. Fit parameters are as in Fig. |7[ The 
agreement is very poor, see text. 
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Fig. 7. Fourier transform of the response function as measured 
(dots) and a fit (solid lines) with the Fourier transform of the 
activity corresponding to the time decay according to Eq. (hol) 
for the photo-oxidized film. Fit parameters are the spontaneous 
decay rate p, the trap concentration c, the trapping probabil- 
ity q, the exciton hopping rate d, and an (irrelevant) overall 
amplitude A. The agreement is quite satisfactory. 



real and imaginary parts with a Fourier transform of the 
activity corresponding to a single exponential in the case 
of the solution, Fig. ^, and to a decay following Eq. ( [lO|) 
in the case of the photo-oxidized and non-odixized films, 
respectively. Figs. and g. 

For the solution. Fig. |6| the experimental data can 
very well be described by a single exponential (Lorentzian 
in Fourier space). The polymer chains have very few de- 
fects, apart perhaps from their end-points which may act 
as traps, and have all the same lengths — this results in an 
almost purely exponential decay of the exciton population 
since the lengths of the trap-free regions are monodisperse 
(in the sense of the decoupled-ring approximation) . 



For the oxidized film. Fig. M, the situation changes en- 
tirely. First, when going from solution to film we introduce 
a coupling between the chains which is not present in the 
solution. Second, due to oxidization, many defects acting 
as traps have been introduced on the film, see Ref. |23|(b). 
Accordingly, Eq. (nO) describes the data well, even though 
it has been devised for one-dimensional systems rather 
than films. This may be explained by the large number of 
traps (around 65% of the effective sites according to the 
fit), which prevent the excitons from 'feeling' the three- 
dimensional nature of the film, as Fig. || clearly illustrated 
in the preceding section. Further support for this effec- 
tive confinement of excitons to a single chain and the 
high number of traps is given in Ref. |23[| (b), based on 
an analysis of the emission and absorption photolumines- 
cence spectra. 

There are, of course, errors in the values of the fit pa- 
rameters. The greatest uncertainty comes from the spon- 
taneous decay rate p, which induces very large errorbars 
on all other parameters as well (not shown in the figures). 
If, however, p is kept fixed at the value indicated in Fig. |7| 
and only the remaining parameters are used for fitting, the 
error in e.g. the trap concentration is found to be about 
15%. Nevertheless it should be kept in mind that our most 
important result here is the functional form of our theo- 
retical curve which is able to reproduce the data rather 
than precise values of the parameters. 

The observed photoluminescence decay-time has been 
greatly increased compared to the solution, in spite of the 
presence of traps. At the same time we know that the ma- 
jor decay path for this sample is non-radiative because the 
photoluminescence quantum yield is only a few percent. 
Hence, the long-living photoluminescence must be the sig- 
nature of a 'stabilized species' which is nof photogenerated 
with a high yield. This may be due to low-energy sites 
on the polymers which capture the excitons and prevent 
them from decaying spontaneously, or delayed photolumi- 
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nescence due to triplet-triplet annihilation. Evidence for 
the latter process in this particular material can be found 
in Ref . [gj| . Neither of these processes have been included 
in our simplified model. 

Turning to the last of the three sets of data, the non- 
oxidized film shown in Fig. pi we observe that it can nei- 
ther be described by a single exponential nor with Eq. (pXJ). 
Fig. M shows the result of an unsuccessful fitting attempt 
with Eq. ([lOl). Apparently, the influence of the three-dimen- 
sional nature of the film made from the same molecules as 
were measured in solution first, cannot be described sim- 
ply by the introduction of traps. Accordingly, since our 
model is not appropriate for the non-oxidized film, the ac- 
tual values of the fit parameters as shown in Fig. || are 
entirely meaningless. 

This has an important implication. The non-oxidized 
film has a photoluminescence quantum yield of 30% vs. 
nearly 100% for the solution of the same molecules. One 
way to account for this could be the formation of traps 
during the film forming process, e.g. by conformational 
stress on the molecules [g5| . If this was the case, our model 
would capture this effect and describe it properly. Hence, 
the reduction of the quantum yield and the difference in 
photoluminescence dynamics between film and solution 
have to originate from a different solid state effect. 



6 Summary and conclusion 

We have considered a model of random walkers undergo- 
ing decay through both capture by imperfect traps and 
spontaneous decay. In one dimension, we have derived an 
explicit analytical form for the time dependence of the 
survival probability, based on the ring-decoupling approx- 
imation scheme, supplemented with an asymptotic expan- 
sion. We have demonstrated that this result, although ex- 
act only in the limit of perfect traps, performs remarkably 
well upon lowering of the trapping probability when con- 
fronted to numerical simulations. Extension of our simu- 
lations to three dimensions, albeit with anisotropic diffu- 
sion, have illustrated the high sensitivity of the decay rate 
to the low-dimensionality constraint. 

An application of our findings to the understanding of 
the decay dynamics of excitons in semiconducting poly- 
mers has been attempted. We conclude that within our 
model traps alone cannot account for the observed dif- 
ference between exciton dynamics in a pristine solution 
and film made from the same conjugated polymer, LPPP. 
Photo-oxidizing LPPP leads to isolated segments which 
probably correspond to the one-dimensional case in our 
model. The large concentration of effective trap sites ob- 
viously dominates over all other physically relevant mech- 
anisms. 

In this context, it is important to emphasize again the 
crucial role played by the spatial trap distribution. For a 
more or less regular spacing of defects, the long-time ki- 
netics would be governed by a simple exponential, with 
a 'renormalized' decay rate. For a random trap distribu- 
tion, on the other hand, a much slower stretched exponen- 
tial decay ensues asymptotically. At least in the long-time 



limit, one might therefore quite drastically control the ex- 
citon population through appropriate 'engineering' of the 
spatial arrangement of the trapping defects. 

At any rate, our model system demonstrates again the 
importance of including statistical ffuctuations effects ef- 
fectively low-dimensional samples. This is of course a well- 
established fact in statistical mechanics model systems. 
Yet to date there have been few examples where clear ev- 
idence of correlations not captured by mean-field theory 
has been found in real experiments on non-equilibrium 
systems. Moreover, we believe that our study underscores 
the relevance of simplified kinetic models, at least in the 
long-time limit, where detailed microscopic mechanism be- 
come less prominent. We hope that a better understanding 
of, e.g., the important organic semiconducting materials 
will eventually be achieved through combining experimen- 
tal data with both quantum-mechanical computations and 
more macroscopic, statistical approaches. 
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A Appendix: Decoupled-ring approximation 



A.l Computation of eigenvectors and eigenvalues 



In this appendix, we show how to compute the symmetric 
eigenvectors and eigenvalues of the matrix 
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ioy 



(12) 



appearing in Eq. (nl). The anti-symmetric eigenvectors, 
also present in principle, are irrelevant in this context be- 
cause they have vanishing overlap with the initial homoge- 
neous particle distribution. We therefore try a symmetric 
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ansatz of the form 



a2 — a 

as 
04 



a„_2 = 04 
\ a„ = 02 = ( 



(13) 



and note that this imphes (Ma)i = a which identifies a 
as the prospective eigenvalue corresponding to a. We can 
then read off that the second component of Ma is given 

by 



(Ma)2 



1 



1 



-as 



(14) 



which must be equal to aa2 since a is supposed to be an 
eigenvector. This gives 



a3 = 2a^ -l+q 



(15) 



By a similar argument, the following recursion relation 
can be derived for the other entries of the eigenvector: 



Uk = 2aak-i - ak-2 (k > 4) 



(16) 



This kind of second-order recursion relation has two expo- 



nential solutions r^ (the —2 in the exponent is for later 
convenience), and ak is given by a linear superposition of 
them, ak = Ar'^~'^ +Br'^~^. Plugging the ansatz r*^'"^ into 
the recursion relation Eq. ( |lq ) yields a quadratic equation 
for r with the two solutions r± = a zt i\J\ — o? . Noting 
that < a < 1 is implied by the non-conserved number of 
particles owing to the trap and the fact that the particle 
density cannot be negative, this can be written in terms 
of an angle (\) with a = cos as rj- = e 



±ii 



The initial conditions 02 = a and 03 = 2a. — 1 + (/ 
now determine the constants A and _B; a straightforward 
calculation gives 



A = B* = ie'* + --^- 
2 e"^ - e- 



(17) 



resulting in 

ai = 1 , 

,, ,,, sin(fc-2),^ ,„ , , (18) 

afc = cosfc-lU + g — ^^ -^— (2<k<n). 

smff) 



= is technically a solution to Eq. dl9), it is not a 
valid one since the constant A does not exist for (f) — 0]. 
Assuming that the smallest </> can be represented by a 
perturbative expansion as = — + ■% + •••, Eq. (|l9| ) can 
be expanded in powers of 1/n and the coefficients Ci are 
readily evaluated. The result is 

(^ - n)'^ 
-Att ^ „ V +Oin-^) . (20) 



^-2.i^ 
n qn^ 



q2j^3 



Overlap of the eigenvectors with the initial state 

As seen in Sec. 13, we also need to know the overlap of the 
eigenvector found in the preceding subsection with the 
initial homogeneous particle distribution, i.e., we need to 
calculate 



(a, 1) = ^flfc , and (a, a) 



fe=i 



k=l 



(21) 



The symbol 1 here stands for the vector with all entries 
equal to 1. Both sums can be evaluated tediously but 
straightforwardly by inserting Eq. (|l8|) and noting that 
both of them can be written in terms of geometric se- 
ries. The resulting expression may be evaluated further 
for large n by inserting cj) = '^/n + 0{n~^). The results 
read 

3 2gn2 



(a,l)" 

71 - 

Altogether this gives 



and (a,a)="^°°^. (22) 



2tt^ 



(a,l)^ n- 
(a,a) 



8n 



(23) 



A. 2 Asymptotic evaluation of the long-time behavior 

Eq. (|l^), the main result of this calculation, can be evalu- 
ated asymptotically for long times. Since (j) = 0{l/n), the 
long-time behavior is dominated by the large-n terms in 
the sum. Using a = cos{<l)) = l — TT'^/2n^+2TT^{l~q)/qn'^ + 
Ojn^'^) and ignoring irrelevant prefactors for the moment, 
Eq. ( fiv( i can be written as 



p{t) ex y^ n exp 

n=l 



"Ml-c)-^i 



2^2(1 -g) 



qn'' 



(24) 



By introducing x = t^'^, the sum on the r.h.s. can be 
recast as a Riemannian sum approximating an integral in 
the variable z = n/x, 



It now only remains to be checked that a„ = 02, as im- 
posed from the beginning. This condition determines p(t) c>c x 
and, after some algebra, results in the following equation: 



^E 



I n ^^ 27r^(l-g) 



_^x[n ln(l — c)/j;- 



X X 



(25) 



[1 — (1 — 2q) cos 20] sin 710 

-[1 - (l-2g) cos 710] sin 20 = . (19) 

The smallest non-zero solution to this equation leads to 
the sought-for largest eigenvalue a = cos [even though 



X — ^CXD 2 

^ X 



dzze2-'(i-'?)/^'«e^[^i"(i-^)-'/2-'l(26) 



This integral may now be evaluated asymptotically by 
Laplace's method and finally yields, including all prefac- 
tors, Eq. (|l]). 
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